Birefringence in nonlinear anisotropic dielectric media 
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Light propagation is investigated in the context of local anisotropic nonlinear dielectric media 
at rest with the dielectric coefficients e''^ = £^^(E,B) and constant ^, in the limit of geometrical 
optics. Birefringence was examined and the general conditions for its occurrence were presented. 
A toy model is exhibited, in which uniaxial birefringent media with nonlinear dielectric properties 
could be driven by external fields in such way that birefringence may be artificially controlled. The 
effective geometry interpretation is also addressed. 

PACS numbers: 04.20.Cv,04.20.-q,11.10.-z,42.25.Lc 



I. INTRODUCTION 

In the regime of intense electromagnetic fields or inside 
some material media the equations governing electrody- 
namic phenomena are nonhnear. In the first case, the 
theory is built from a nonlinear Lagrangian which is a 
function of the two Lorentz invariants of the electromag- 
netic field. In the second case. Maxwell equations must 
be supplemented with constitutive relations between ex- 
ternal applied fields and the induced excitations. In gen- 
eral, such relations are nonlinear (although linear consti- 
tutive relations are also of interest "5] ) , and they depend 
on the physical properties of each considered medium un- 
der the action of external fields. In both cases the field 
equations will be presented in a nonlinear form and, as a 
consequence, several non usual effects (in the context of 
Maxwell theory) are predicted. 

Of particular interest is the phenomenon of birefrin- 
gence: the light velocity dependence on the polarization 
mode of the propagating wave. In the case of nonlin- 
ear electrodynamics, birefringence occurs when the elec- 
tric and/or magnetic fields exceed its critical value pre- 
dicted by quantum electrodynamics The analysis of 
light propagation in such regime shows that there is a 
probability of photon splitting under a strong external 
electromagnetic field More recent investigations on 
light propagation in the context of nonlinear Lagrangian 
for electrodynamics can be found in Refs. IE M H Hi 
and references therein. Particularly, in Ref. 8] the Fres- 
nel analysis of wave propagation in nonlinear electrody- 
namics was performed. For a class of local nonlinear 
Lagrangian nondispersive models, birefringence phenom- 
ena were studied and the wave propagation in a moving 
isotropic nonlinear medium was also examined. 

Natural uniaxial birefringence is a well known effect 
and it takes place for some materials (mostly crys- 
tals). Artificial birefringence is also possible to occur 
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[Tfil ITU IT^ , as an induced effect in material media: 
when an external field is applied in a medium with non- 
linear dielectrics properties, an artificial optical axis may 
appear. Birefringence is also used as a technique for in- 
vestigating other properties of some systems: see for in- 
stance Ref. where the birefringence effect is used as 
a tool for astrophysical studies. 

In this work, birefringence is investigated in the con- 
text of homogeneous dielectric media at rest with the 
dielectric coefficients e^i, — e^,y{E,B) and constant fi 
in the limit of geometrical optics. The analysis is re- 
stricted to local electrodynamical models, and disper- 
sive effects were neglected by considering only monochro- 
matic waves, thus avoiding ambiguities with the velocity 
of the wave. A mechanism by which birefringence can be 
controlled by means of an external electric field is pro- 
posed. In particular, it is shown that naturally uniaxial 
media presenting nonlinear dielectric properties can be 
operated by external fields in such way that birefringence 
could be artificially turned off. 

In the following section, the generalized eigen-vector 
equation associated with the light propagation in general 
anisotropic media is presented. In Sec. lIIII the solution for 
such equation is presented for the particular case where 
the impermeability fx^i, is assumed to be a constant, and 
birefringence conditions are examined. Section llVl partic- 
ularizes to the case e°'p{E), and the principal refractive 
indices are explicitly presented. In Sec. a particular 
model is exhibited in which a naturally uniaxial media 
presenting nonlinear dielectric properties, under the ac- 
tion of an external electric field, behaves in such way that 
birefringence could to be artificially controlled. 

A covariant formalism is used throughout this work. 
Spacetime is assumed to be Minkowskian, and a Carte- 
sian coordinate system is used, such that the metric is 
'7^1/ — diag(+l, —1, —1,-1). Units are chosen such that 
c = 1, except mentioned otherwise. A geodetic observer 
V^^ = 6q is supposed to describe all quantities. Particu- 
larly the electric field is represented by E^^ = —F'^'^Vi, ~ 
(0, E) whose modulus is E ~ {—E^EaY^'^ , and similarly 
for the magnetic field B^' = (0, B). 
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II. WAVE EQUATION 

The electrodynamics in a medium at rest is completely 
determined by the Maxwell equations 



(1) 
(2) 



together with the constitutive relations 

L>" = e'^fjE'^ (3) 
i/" = A^^^B^ (4) 

where the coefficients e"p = e"p{E, B) and /x"/3 = 
li°'f}{E, B) represent the dielectric tensors (called per- 
mittivity and impermeability tensors, respectively) which 
encompass all information about the dielectric proper- 
ties of the medium. They are usually assumed as being 
frequency-dependent, although they may more generally 
depend on the external fields themselves. The coordinate 
derivatives of Eqs. Q and Q are 



dE- 



^E^^B^, (5) 



^ dB^ 



dE^ 



B^E 



In order to determine the propagation of the electro- 
magnetic waves, we will consider the eikonal approxima- 
tion of electrodynamics, making use of the method of field 
discontinuities jlj, il5> | . With the notation introduced in 
Ref. nil, we set [E^' ^^]^ = e^'K^ and [B^' ^y]-^ = b^'K^, 
where e'^ and 6^ represent the polarization of the wave 
fields on the surface of discontinuity S and K\ = d'S/dx'^ 
is the wave vector normal to S. By applying these bound- 
ary conditions to the Eqs. ^ and Q, and taking into ac- 
count Eqs. (0 and 10) , we obtain 6'^ = if^^K^Vxe^/u, 
and 



VxX.rV^K^V'Kp = 0, (7) 



where we introduced the new quantities 



(9) 



with uj = K^Va the frequency of the electromagnetic 
wave. 

In order to present Eq. Q in a 3-dimensional rep- 
resentation we introduce the projector on the 3-space 



h'^ = S'i^- V^V^, and define 



(?" = h'^Rt' = K°' - ujV° 



(10) 



Straightforward simplifications lead Eq. ^ to the form 
Z°'re^ = 0, (11) 

where 

Z\ = C\ + l^rj'^f^-'^B'V.qf, + ^I^;I^^HA12) 

in which a^^'^l = a^'^ — a"^ for an arbitrary tensor a^^ . 
Further, we have defined the projector 



(13) 



on the 2-space orthogonal to the propagation direction 
qt^ J- qt^ /q^ where = —h^uq^q"- 

The general solution for the wave propagation can be 
derived from the eigen- vector problem stated by Eq. l(TT|l. 
and is given by the generalized Fresnel equation 



det|Z"^| = 



(14) 



For several physical configurations 6] the dispersion re- 
lations coming from Eq. (|14|l can be written in the sug- 
gestive form g'^K^K^ = 0, as shown in Sec. IIIII 



III. BIREFRINGENCE CONDITIONS 



Let us assume the more specific case where the medium 
behaves in such way that /x"^ = ^~^hp with ^ a constant. 



In this case, Eq. p2l) reduces to 



(15) 



where ~ uj/q is the phase velocity of the wave. Equa- 
tion ini) can be written as (Z''^)^ - SZ^^Z"^^''^ + 
2Z°'i3Z^~^Z^a = 0, which yields a fourth order equation 
of the form 



av^ - /3w<^ - 7 = 0, 
where we have defined 



(16) 



" = J [(C%)3-3C%C"0C^+2C%C^C^„],(17) 



P = f^-' {C^cC^" - C^C^") qxq., 
7 - fi-^C^-'qxq,. 

The solutions of Eq. (|16|l are 




1± A 1 



4a7\ 



(18) 
(19) 



(20) 



Equation (|20|) expresses the fact that, in general, the 
velocity of the electromagnetic waves inside a material 
medium may get two possible values w J , v~ which are 
associated to the two possible polarizations modes T^. 
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Thus, birefringence occurs if the following inequalities 
hold for a given direction q'^: 



aP > 0. 



(21) 
(22) 



(For the particular case e'^p — eft-^, with e a constant, 
one obtains a = e^, (3 = 2e^/^ and 7 = — e/^^, and the 
two phase velocities reduce both to the same well-known 
value v"^ = l/efj,.) 

A. Effective geometry 

The dispersion relation Eq. (|16|l can be conveniently 
expressed in the form gj^K^K^ ~ (corresponding to 
the factorization of the general fourth-order Fresnel ten- 
sor ^3), where 



a/3 

9± 



(23) 

The symmetric tensors g'^ represent the effective geome- 
tries (also known as optical metrics) associated with the 
wave propagation, and the symbol ± indicates that, in 
general, there will be two possible distinct metrics, one 
for each polarization mode. The integral curves of the 
vectors are geodesies in the corresponding effective 
geometry gj^ 12]. For the particular case of the linear 
theory in vacuum, both g'^ and g'^ reduce to ry^" as 
expected. For the case of naturally isotropic nonlinear 
media at rest the results presented in Refs. (3. lllllT^ are 
recovered. 

The effective geometry interpretation for electromag- 
netic waves in material media 0, 0, 0| (or even 
in the context of nonlinear electrodynamics in vacuum 
lEIE^I) has been proposed as a tool for testing kine- 
matic aspects of general relativity in laboratory '19^. In 
this context, Eq. (|23|1 provides a natural support for 
analogue anisotropic models. 



IV. NON-MAGNETIC ANISOTROPIC 
NONLINEAR MEDIA 

Let us consider a naturally uniaxial medium that reacts 
nonlinearly when subjected to an external electric field as 
e"^ = diag[0,£||(£'),e_L(i?),£^(i?)]. In this case, e"^ — 
e" 13(E) and Eq. (jSJ takes the simpler form 



^ 13 - £ t 



~dE 'e~ 



(24) 



By setting E in the x-direction (optical axis), we obtain 
C"' p = diag(0,e|| + i?£'|,ej_,£_L), where £|| = de\\/dE. 
Now, from the above results, we obtain 



07 



(25) 



4a7 -f- fr 



(26) 



C\ [C\{1 
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£i(l-<7i')] >0, (27) 



which show that conditions (|21|l and (|22|l are fulfilled in 
this case, provided that C^i 7^ £± and qi^ < 1, since 
both C^i and e± are positive quantities. 

Equation H26|l shows that birefringence can be annulled 
if C^i = £_L. In other words, birefringence phenomena 
could be turned off by imposing that 



■Ee[ 



(28) 



which is a no-birefringence condition. Note that if 
a model for dielectric permittivity is set, the no- 
birefringence condition will mean a condition on the value 
of the external electric field. 



A. Ordinary and extraordinary rays 

By considering the previous model, Eq. H2U|) yields 
1 



2ne±Ch 
±{C\-e^)il-q^)] 



(29) 



Thus, there will be two phase velocities, one of which is 
isotropic (ordinary ray) 

1 



/i£j 



(30) 



and the other one depends on the direction of propaga- 
tion (extraordinary ray) 



1 



[e^{l-qi^) + C\qi' 



(31) 



The two velocities coincide when either the propagation 
occurs along the direction of the electric field, or when 
the no-birefringence condition Eq. (|28|l holds. 

The two extreme values of v~ occur in the direction 
of the electric field {9 = 0), and perpendicularly to it 
{9 = tt/2). For these two directions, we obtain 



^11 



le=o ^£^ ^ ^ 
1 

= (V)'L=./9 = 



(32) 
(33) 



Associated to these directions we define the principal ef- 
fective refractive indexes and n^^ by 



1 



/i£j 



(34) 
(35) 



The effective refractive index of the extraordinary ray 
along an arbitrary direction q^ is determined by Uq — 
with v~ given by Eq. (|^ . 
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CONTROLLED BIREFRINGENCE 



Let us examine some special cases where birefringence 
could be controlled by adjusting the magnitude of an 
external electric field. 



A. Non-birefringent nonlinear media 

We consider the model where = — 3pE^, with 
e_L and p constants. Thus, from Eq. (|28|l we obtain 

equation at i? = yields 



-Ee'l — e±~ 3pE^. Regularity on the solution of this 



pE^ 



(36) 



It means that for every value of E, if the medium reacts 
like Eq. there won't be birefringence phenomena. 



B. Anisotropic media with artificially controlled 
birefringence 

Now, let us examine the model for which the medium 
reacts as 



= ej_-3pE^ 
= e|| - sE^, 



(37) 
(38) 



with ej^, p, e|| and s constants. For the case p = s, all 
values of E will correspond to either a no-birefringent 
configuration (if e± = e||), as described by Eq. H3()|) . or 
else a bircfringcnt configuration (if e± ^ e||). On the 
other hand, for p ^ s, Eq. (^5)1 yields e|| — 3sE^ — 
e± — 3pE^, and solving it for the electric field one obtains 



Er_ 



3(s-p)- 



(39) 



Thus, from Eq. H39|) we conclude that there will be a 
particular value Ec of the electric field for which the no- 
birefringence condition will be satisfied. Any other value 
of E Ec will result in birefringence. Of course these 
results are only of interest in the case Ec lies within the 
limit of applicability of Maxwell theory. Further, the 
quadratic terms in E in Eqs. (|37|l and H38(l should remain 
small when compared with the constant ones. Equation 
is meaningful for positive media {£± < if p < 
and s > 0; and similarly for negative media (e_L > if 
p > and s < 0. 



For the case where ey ^ e± birefringence will occur pro- 
vided that E Ec- For the particular case where E — 
we have natural birefringence, as it occurs in some crys- 
tals. See Fig. ^for a numerical analysis of a toy model 
where natural birefringence is artificially controlled by 
the influence of an external electric field. 




FIG. 1: Controlled birefringence: the particular value Ec of 
the electric field for which birefringence is turned off occurs 
in the intersection of the curves representing the effective re- 
fractive indices for the ordinary and extraordinary rays. The 
naturally uniaxial positive crystal becomes negative when the 
electric field exceeds Ec- Note that, for E = 0, natural bire- 
fringence remains. We used the following illustrative values 
(in MKS units) for the parameters appearing in Eqs. 13711 
and lEU: ey = 1.5195 x 10""F/m, £± = 1.5172 x 10""F/m, 
p = -2.2648 X 10~"Fm/V^ and s = 2.2700 x lO^^'^Fm/l/^ 



By considering a naturally isotropic medium where 
e|| = e_L, the model set by Eqs. (jTzjl and describes ar- 
tificial birefringence ^ it appears in the Kerr electro- 
optic effect. For this case we obtain Uo — ric w XKE^, 
where we have defined XK = 3y/JIe±{s —p)/2e±, with 
A the wave length of the electromagnetic waves and K 
the Kerr's constant (its value depends on the dielectric 
properties of the medium). The effective geometry for 
this particular case has already been presented in the 
literature HEiilJ. 



VI. CONCLUSION 

A covariant tcnsorial formalism is here provided in or- 
der to discuss the propagation of monochromatic elec- 
tromagnetic waves inside naturally anisotropic material 
media with nonlinear dielectric properties, in the limit of 
geometrical optics. The eigen- vector problem for general 
media was presented. For the case of constant /i, the con- 
ditions for birefringence phenomena to appear were ex- 
amined. For the case of non-magnetic media, a particular 
model was shown for which the magnitude of the bire- 
fringence could be controlled, and even turned off, with 
the use of an external electric field. The model general- 
izes the standard description of the Kerr electro-optical 
effect to the case of a naturally anisotropic material me- 
dia. All the obtained results can be straightforwardly 
rewritten to the case e"p = ehp with e = const and 

M"/3 = M"/3(^); thus generalizing the Cotton-Mouton ef- 
fect (magnetic analogue to the Kerr effect) to the case of 
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anisotropic media. 

The results can also be used in the context of analogue 
models for general relativity. It is here presented the 
optical metric for nondispersive anisotropic media, which 
extends previous investigations. Isotropic moving media 
were already discussed and their results agree with 
ours in the case of isotropic media at rest. Anisotropic 
moving media are still to be considered. 
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